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ABSTRACT 

We have used fibre spectroscopy to establish cluster membership and examine pre- 
main-sequence (PMS) lithium depletion for low-mass stars (spectral types F to M) 
in the sparse young (2± 30Myr) cluster IC 4665. We present a filtered candidate list 
of 40 stars that should contain 75 per cent of single cluster members with 11.5 < 
y < 18 in the central square degree of the cluster. Whilst F- and G-type stars in 
IC 4665 have depleted little or no lithium, the K- and early M-type stars have depleted 
more Li than expected when compared with similar stars in other clusters of known 
age. An empirical age estimate based on Li-depletion among the late-type stars of 
IC 4665 would suggest it is older than 100 Myr. This disagrees entirely with ages 
determined either from the nuclear turn-off, from isochronal matches to low-mass 
stars or from the re-appearance of lithium previously found in much lower mass stars 
(the "lithium depletion boundary"). We suggest that other parameters besides age, 
perhaps composition or rotation, are very influential in determining the degree of 
PMS Li-dcplction in stars with M > 0.5Mq. Further work is required to identify and 
assess the effects of these additional parameters, particularly to probe conditions at 
the interface between the sub-photospheric convection zone and developing radiative 
core. Until then, PMS Li depletion in F- to early M-type stars cannot be confidently 
used as a precise age indicator in young clusters, kinematic groups or individual field 
stars. 
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1 INTRODUCTION 

The lithium depletion boundary (LDB) technique deter- 
mines ages for young clusters of stars by establishing the 
luminosity at which Li remains unburned in the atmospheres 
of fuUy-convective, very low-mass stars. In principle, LDB 
ages are both precise and accurate; observational and the- 
oretical uncertainties contribute to errors of no more than 
about 10 per cent in the age estimate for clusters in the 
range 10-200 Myr (Jeffries & Naylor 2001; Burke et al. 2004) 
- considerably better than other age estimation methods. 

LDB ages have now been estimated for several clus- 
ters (Stauffer et al. 1998; Barrado y Navascues, Stauffer & 
Jayawardhana 2004; Jeffries & Oliveira 2005) , but the tech- 
nique has limited applicability because it entails quantify- 
ing the Li I 6708A feature at a resolving power R ^ 3000 
in groups of faint, very low-mass objects. However, the few 
LDB ages that are known could be used to calibrate other 
age estimation methods that are feasible in more distant 



populations and isolated field stars, but which rely on more 
uncertain stellar physics. 

A prime example is the age-dependent degree of Li 
depletion seen in more massive (and many magnitudes 
brighter) F-K stars. Li is burned in pre-main-sequence 
(PMS) stars as their contracting cores reach temperatures 
^ 2.5 X 10^ K. In fuUy-convective stars this leads to rapid and 
almost total Li destruction on less than a Kelvin-Helmholtz 
timescale. However, in the standard picture, PMS stars with 
M > 0.5 Mq form radiative cores that push outwards, so 
that the base of the convection zone dips below the Li igni- 
tion temperature, bringing an abrupt halt to surface Li de- 
pletion (e.g. Piau & Turck-Chieze 2002). The exact details 
and in particular, isochrones of Li depletion versus mass or 
effective temperature, are extremely sensitive to conditions 
at the base of the convection zone and any additional mix- 
ing into the radiative core. As a result, whilst Li depletion 
among F-K PMS stars has a long history of use as an empir- 
ical age indicator among young clusters, kinematic groups 
and even individual stars (e.g. Duncan & Jones 1983; Stauf- 
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Figure 1. A V vs y — / color magnitude diagram from our plio- 
tometric survey of IC 4665. Large symbols indicate targets (se- 
lected solely from this diagram) that were spectroscopically ob- 
served at the WHT. A 30Myr isochrone from the Baraffe et al. 
(2002) models is shown (using a colour/effective temperature re- 
lationship that is tuned using data from the Pleiades - see Jeffries 
& Oliveira [2005] for details), adjusted for a distance of 370 pc, 
Av = 0.59 mag and E{V - I) = 0.23 mag. 

fer, Hartmann & Barrado y Navascues 1995; Favata et al. 
1998; Song et al. 2000), the absolute ages determined in this 
way have large theoretical uncertainties. 

In this paper we discuss results from a programme to 
calibrate Li depletion ages among F-K stars by observing 
such stars in clusters with known LDB ages. We present pho- 
tometry and spectroscopy of F-K (and M) stars in IC 4665, 
an open cluster with a distance of 370 ± 50 pc, a reddening 
E{B - I/) = 0.18 ± 0.05 mag (Hogg & Kron 1955; Crawford 
& Barnes 1972 - corresponding to E{V - I) = 0.23 ± 0.06 
and Av = 0.59 ±0.16, for the intrinsic colours of a low-mass 
PMS star [Bessell et al. 1998]) and an unambiguous detec- 
tion of the LDB which yields an age of 28 ± 5 Myr (Manzi et 
al. 2008, hereafter M08), where the modest age uncertainty 
is dominated by the distance uncertainty. 



2 OBSERVATIONS AND ANALYSIS 

Spectroscopic observations of candidate low-mass stars in 
IC 4665, with 11.5 < F < 18 and 0.4 < V - / < 3, were 
conducted on the nights beginning 8 June 2008 and 9 June 
2008, using the AF2/WYFFOS muhi-fibre spectrograph at 
a Nasmyth focus of the 4.2-m William Herschel Telescope 
(WHT). An echelle grating and an order-sorting filter were 
used to give a spectral coverage of about 420A centred at 
6600 A and a resolving power of 11 000. 

Targets were selected from a ~ 1 square degree photo- 
metric survey centred at RA= 266.55°, Dec=-|-5.7° (J2000), 
performed on 15-17 September 2005 using the 1.0-m tele- 
scope at the Cerro Tololo Interamerican Observatory. The 
details of this survey will appear elsewhere (Cargile et al. 
in prep.), but in brief consisted of nine 19.3 x 19.3 arcmin^ 
fields observed for 100, 50 and 25 seconds through B,V,I 
filters. The observations were reduced with daophot (Stet- 
son 1987) and placed onto the Johnson- Cousins system us- 
ing standards from Landolt (1992). The V,V — I colour- 
magnitude diagram (CMD) for the survey is shown in Fig. 1, 



along with targets observed at the WHT. Targets were se- 
lected from this diagram assuming the LDB age of 28 Myr, 
but the criteria were broad enough to encompass all likely 
cluster members and all previously identified low-mass mem- 
bers in this magnitude range (e.g. Prosser & Giampapa 
1994 - hereafter PG94; Martm & Montes 1997 - hereafter 
MM97). A total of 452 targets were observed in 8 fibre con- 
figurations, amounting to 87 per cent of possible candidates 
with 11.5 < < 18. Four "bright" configurations were ob- 
served for 3 X 600 s and 4 "faint" configurations were ob- 
served for 3 X 1800s (one had only two good exposures). 
The spectra had signal-to-noise ratios (SNR) ranging from 

4- 110 per O.22A pixel. 124 targets were observed twice, 37 
targets were observed three times and 6 targets were ob- 
served in four configurations. 

Spectra were extracted from the CCD images using pur- 
pose built software (written by RJJ) which used an optimal 
extraction procedure to maximize the signal to noise ratio. 
Science images were corrected for bias by subtracting a me- 
dian of multiple bias frames and correcting for any time- 
dependent offset using overscan regions. Bad pixels were 
flagged using a mask generated from dark frames. Tungsten 
lamp images recorded before each science run were used to 
trace the path of the spectrum on the CCD image. The me- 
dian of these images was used to define flat-field spectra for 
the optimal extraction process. Copper-neon lamp spectra, 
recorded before and after science exposures, provided the 
wavelength calibration. 

The method used for optimal extraction is described 
by Horne (1986). This applies non-uniform weight to pix- 
els in the extracted sum, minimizing statistical noise whist 
preserving photometric accuracy. In our case the weighting 
profiles for individual spectra were derived from a boxcar 
average along the wavelength axis of the flat-field images. 
A similar boxcar averaging process using data from regions 
between the spectra was used to model and subtract any 
scattered light in each exposure. With this procedure the 
extracted science spectra were normalized to the flat-fleld 
spectra, compensating for variations in gain between pix- 
els. To reduce the effect of cosmic rays the measured signal 
profile was compared with the product of the extracted spec- 
trum and the weighting proflle. Pixels showing greater than 

5- sigma difference were masked and the extraction process 
repeated. Profiles with more than 50 per cent of the pixels 
masked were masked in the flnal extracted spectrum. 

Arc spectra, extracted using the same procedure, were 
used to identify the position of seven well separated arc lines 
and a cubic polynomial flt was used to deflne the calibration 
of individual spectra. Extracted spectra and their variances 
were interpolated across masked points and recast onto a 
common wavelength base comprising 10 000 logarithmically 
spaced steps between 6350A, and 6740A. A median sky spec- 
trum, determined from ten or more flbers pointed at sky po- 
sitions in each configurations was subtracted from the target 
spectra within that conflguration. Spectra of the "blank" 
night sky were taken to measure the relative efficiency of in- 
dividual fibers in each conflguration. The results of repeated 
exposures were averaged to produce the median spectrum of 
each target from each conflguration. 

Spectra of late-type stars with known radial velocities 
(RVs) (from Nidever et al 2002) were taken on each night. 
Heliocentric RVs were measured by cross-correlation using 
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the IRAF FXCOR task and HD 171067 (spectral type G8V) 
as the principal template. The cross-correlation used the 
wavelength regions 6380-6550A and 6580-6725A. Projected 
equatorial velocities (usini) were estimated from the width 
of the cross-correlation function. Translation from the cross- 
correlation function width to a w sin i was calibrated by ar- 
tificially broadening standard star spectra. Variation in the 
intrinsic resolution of difi'erent fibres was accounted for us- 
ing the cross-correlation function between IfD 171067 and 
the twilight sky spectrum recorded through each fibre. At 
this resolution the data were capable of resolving rotational 
broadening corresponding to v sin i > 20 km s~ ^ . FXCOR pro- 
vides error estimates for the RVs, which concurred with 
estimates from targets observed in multiple configurations. 
The average RV uncertainty was 0.8 km s^^, but varied from 
0.3kms~^ to as large as 8kms^^ in the most rapid rotators. 
RVs were obtained for 449/452 targets, the remaining 3 tar- 
gets had low SNR and no discernible cross-correlation peak. 
V sin i errors were estimated by numerical simulation and 
were of order 10 per cent, where rotational broadening was 
resolved. 

The equivalent widths of the Li I 6707. 8A feature 
(EW[Li]) were measured with respect to a pseudo- 
continuum, using a Gaussian fit to the line itself where 
present. Uncertainties were estimated using the empiri- 
cal SNR of the continuum fit and the formula crEW = 
1.5y^pw/SNR, where p is the pixel size and w the FWHM of 
the line in wavelength units (Cayrel 1988) . Where no signif- 
icant line was seen we estimated very conservative 3-sigma 
upper limits using this formula and an assumed w (including 
any rotational broadening). At our resolution, the measured 
EW[Li] includes a blended contribution from a weak (~ 10 
to 20mA) Pel line at 6707.44A. 

A summary of the results is included as Table 1 (avail- 
able fully in electronic form only). This identifies each target, 
gives its J2000 celestial coordinates (taken from the 2MASS 
catalogue - Cutri et al. 2003), its V and V — I magnitudes, 
the heliocentric RV measurements, a weighted mean RV, a 
RV variability fiag, a flag identifying objects that are candi- 
date SB2 binary systems based on the appearance of the 
cross-correlation functions, the iisini (or upper limit), a 
SNR for the (summed) spectrum and a EW[Li] (or upper 
limit) found from the (summed) spectrum of the object. Ta- 
ble 1 also lists alternative star designations from Prosser 
(1993) and M08. Figure [2] shows a sample of our spectra 
around the Ha and Li I 6708A features. 



3 CLUSTER MEMBERSHIP 

Cluster membership was assessed primarily from the RVs. 
The cluster is clearly present as a concentration in the 
RV distribution (see Fig. [Sjl, although it sits upon a 
much broader contaminating field star distribution. The 
subset of stars with strong Li 6708A features have very 
similar RVs. Discarding several outliers, there are 22 Li- 
rich stars, which we assume are cluster members, with 

— 15.2 <RV< — lO.lkms"^, giving a weighted mean RV of 

— 12.9 ± 0.3kms~^. Using spectra of all the observed RV 
standards we estimate an additional external error on this 
mean of ±0.5kms~^, making our value consistent with the 
-12.9 ± 0.4kms-i found by PG94, but 3.0 ± 1.2 kms'^ 
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Figure 3. The measured equivalent width of the Li I 6708A fea- 
ture aginst heliocentric radial velocity. Solid points are detections, 
whereas triangles represent 3-sigma upper limits. There is a clear 
concentration of Li-detections at an RV ~ — 13kms~^. 
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Figure 4. The summed equivalent width of the Cal lines at 
6349A and 6463A as a function of colour. Candidates with a de- 
tected Li I feature or Ho in emission are indicated. The candidates 
in the box have weak Cal and we class them as probable giants 
and hence non-members. 



higher than the value reported by M08. Comparing the stan- 
dard deviation of the RVs with their individual uncertain- 
ties indicates an intrinsic cluster RV dispersion of 1 kms^^. 
Adding this internal dispersion to the RV uncertainties in 
quadrature, we then accept any star with an RV within 1.5- 
sigma of the cluster mean as a probable cluster member, 
yielding 56 candidates. Applying similar criteria for assumed 
mean RVs 10 km s~^ either side of the IC 4665 RV allows us 
to estimate that about 16 ±3 of these 56 candidates are non- 
members with RVs that randomly coincide with IC 4665. 

The candidate list was refined in several ways. Fig- 
ure 1 shows that the heaviest contamination is expected for 
1.0 < V — J < 1.5 where a "finger" of background K-giants 
cuts across the cluster sequence. To filter these out, we mea- 
sured the EWs of two gravity-sensitive, isolated Cal lines 
at 6439A and 6463A. We have used the spectral synthesis 
software uclsyn (Smalley, Smith & Dworetsky 2001) along 
with Atlas stellar atmosphere models (Kurucz 1993) to es- 
tablish that these should be a about a factor of 2 weaker in 
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Figure 2. Examples of spectra in the vicinities of the Ho and Li I 6708A lines. The spectra have been normalised to unity and offset 
for clarity. These examples show the full range of signal-to-noise ratios and spectral types for some candidate members of IC 4665 (see 
Table 2). 



giants than in dwarfs of the same temperature/colour. Ac- 
cordingly, in Fig. 13] we identify 10 of our candidates, all with 
1.2 < V — I < 1.33, that have weaker Cal lines than ex- 
pected for dwarfs and which we class as probable giants (or 
dwarfs with a low metallicity) and hence not cluster mem- 
bers. As expected, none of these have measurable Li, rapid 
rotation or Ha emission (see below). 



Proper motions for 45 of the 56 RV candididates are 
found in the NOMAD database (Zacharias et al. 2004). The 
mean proper-motion of the Li-rich stars is —0.7 ± 1.0 milli- 
arcsec/yr in RA and —6.2 ± 0.8 milli-arcsec/yr in Dec. This 
is not greatly dilferent from the mean for field stars in this 
direction. We find that 2 stars (the brightest F-star in our 
sample [JC04_295] and an M-dwarf with Ha in absorption. 
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Table 1. The spectroscopic measurements of 452 photometric candidates in IC 4665. Only a few rows are shown here to illustrate the 
content, the full table is available in electronic form. Column 1 is an identifier from the photometric catalogue (Cargile et al. in prep), 
columns 2 and 3 give the RA and Dec (J2000, from the 2MASS catalogue), columns 4 and 5 give the V and V — I photometry, column 
6 lists how many observations were made of the target then columns 7-9, 10-12, 13-15 and 16-18 give up to 4 possible measurements of 
the heliocentric RV, its error and the HJD of the observation. Columns 19 and 20 give the weighted mean RV and its uncertainty, column 
21 gives the mean signal-to-noise ratio of the averaged spectrum, column 22 and 23 give the measured EW[Li] and its uncertainty. An 
uncertainty of -9999 indicates that the quoted value is a 3-sigma upper limit. Columns 24 and 25 list the projected equatorial velocity 
and its uncertainty. An uncertainty of -9999 indicates that the vsini is an upper limit. Column 26 indicates whether the object has 
a variable RV {= 1 if variable at > 99.9 per cent confidence) and column 27 indicates whether the cross-correlation function with the 
template standard star shows evidence for more than one star {= 1 for evidence of binarity). Columns 28 and 29 list alternative names 
for the targets from Prosser (1993) and Manzi et al. (2008). 
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the 4 cool targets {V — I > 2) with Ha in absorption are 
members, because Ha emission is ubiquitous for such stars, 
even in older clusters like IC 2391 and the Pleiades (see 
Stauffer et al. 1997). These four objects are probably older 
contaminating field stars and one of them is also a non- 
member by proper motion as explained above. We -wish to 
use Ha emission as an age indicator later on (see section O, 
so -we do not discard these objects, but mark them as "pos- 
sible" members. 

Finally, -we have taken JHK magnitudes from 2MASS 
(Cutri et al. 2003) and plotted the J — H vs H — K diagram. 
Only one star stands out from the general trend of photo- 
spheric colours, and this only by 2.5-sigma. JC05_329 is a 
target with J- H = 0.15 ± 0.04 and H - K = 0.15 ± 0.04, 
which is too red in H — K hy about 0.10 mag compared 
with other cluster candidates. This could perhaps be due 
to warm circumstellar dust, although this is unlikely at an 
age of ~ 30 Myr. For now this star is classed as a "possible" 
member. 

The 40 stars that pass all these tests are classified as 
definite members. The full 56-star candidate list and a sum- 
mary of the measurements and membership judgements are 
shown in Fig. [5] and listed in Table 2. 



3.1 Comparison with other work on IC 4665 

Many of our targets and members were selected as such and 
studied by Prosser (1993), PG94, MM97 and M08. We com- 
ment in detail on the comparison of our measurements and 
membership selections with these papers below. 
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Figure 5. AV vsV — I color magnitude diagram for the 56 stars 
in IC 4665 with RVs compatible with cluster membership. Stars 
which have Li in absorption or Ho emission are indicated, as are 
12 stars judged to be non-members on the basis of their giant-like 
gravities or discrepant proper-motions (see text). Baraffe ct al. 
(2002) isochrones at 20, 30, 50 and 100 Myr are also shown (these 
use a Pleiades-tuned colour /effective temperature relationship - 
see Jeffries & Oliveira [2005]). 

[JC04_558]) have proper- motions incompatible with cluster 
membership at the 3-sigma level. 

We measured the EW of the Ha line by direct inte- 
gration above (or below) a pseudo-continuum. Most of the 
cooler targets exhibit Ha emission (see Fig. |4]), as expected 
for young, chromospherically active stars. It is unlikely that 
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Table 2. Properties of candidate RV members of IC 4665. Columns 1-6 identify the star and give its photometry, EW[Li] , RV and v sin i. 
An uncertainty of -9999 indicates an upper Umit. Column 7 gives the EW of the Ha feature (positive = emission). Columns 8 and 9 give 
the EWs of the gravity-sensitive Cal lines and columns 10 and 11 give the proper motion from the NOMAD catalogue. Column 12 gives 
our final assessment of membership (see section |3j . 
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We find 137 of our targets have V and / photometry in 
Prosser (1993). The agreement of this photometry with that 
presented here is good. The average difference in V magni- 
tudes (ours minus Prosser) is —0.028 ± 0.007 mag with a 
standard deviation of 0.086 mag. After converting Prosser's 

V — I Kron colours to V — I Cousins colours using the for- 
mulae of Bessell & Weis (1987), the average difference in 
colour is —0.033 ± 0.007 mag with a standard deviation of 
0.077 mag and no obvious colour dependence in the range 
0.5 < - / < 3. 

There are 30 stars with RVs common to our work 
and PG94. The agreement is excellent apart from 4 ob- 
jects. Three of these (JCO4_057=P166 , JCO5_160=P19 and 
JCO8_038=P38) have been identified by us as probable SB2 
binary systems (see Table 1). The other (JC09_462=P66) 
has a RV discrepancy (our value minus that of PG94) of 
— 17.2 ± 0.5kms~ and it is probably a short period SBl 
binary. For the remaining 26 stars the average RV dis- 
crepancy is —0.25 ± 0.41 kms~^ with a standard deviation 
of 2.1kms"\ I gnoring the SB2 objects mentioned above, 
then the v sin i values quoted by PG94 also agree well with 
our measurements. The one exception is JC06_347=P12, 
where PG94 find iisini < lOkms"^ but we have vs'mi = 
108±4kms"^ As both MM97 and AUain et al. (1996) also 
find that P12 is a rapid rotator, with v sin i = 70 km s~^ and 
a rotation period of 0.60 days respectively, then it seems 
likely that the PG94 measurement is in error. We have 12 

V sin i measurements for targets in common with Shen et al. 
(2005) and all are in perfect agreement. 

There are 59 targets common to this paper and M08. 
A comparison of RVs reveals three probable SBl bina- 
ries (JCO5_039 = number 103 in M08, with RV differ- 
ence of 24.2 ± 0.8 km s'^ JC04_271 = M08 number 40, 
ARV= 9.2 ± l.Okms-^ JC04_249 = M08 number 116, 
ARV= 9.9 ± l.Okms"^). After clipping out these objects, 
the mean difference between our RVs and those of M08 is 
2.9 ± 0.3kms~^, with a standard deviation of 1.9kms~^. 

We have been able to compare our values for EW[Li] 
with 12 stars also studied by MM97. Agreement is good 
with the average EW difference (ours minus MM97) being 
— 18 ± 9 mA with a standard deviation of 27 mA. We also 
have 5 targets with EW[Li] listed by M08 and again, the 
agreement is good. 

In terms of selected members, there are 7 objects se- 
lected by PG94 or by MM97 as cluster members that are 
not in our list of candidate members (Table 2). For P12, 
P19, P27, P38 and P166 these either have an RV outside our 
selection range or are possible SB2 binary systems. The ob- 
jects P39 and P155 may be genuine members but are among 
the 13 per cent of photometric candidates we did not observe 
spectroscopically. There are 27 candidate members consid- 
ered by M08 that do not appear in our Table 2. Most of these 
are too faint {V > 18) to have been included in our sample, 
but six (identified by M08 with the numbers 6, 16, 20, 59, 83 
and 125) were observed by us but had RV values outside our 
selection range. JCO5_190 {— M08 number 59) is a proba- 
ble binary system according to both us and M08. JCO6_038 
{— M08 number 20) is a binary according to M08, but we 
have only one RV measurement. M08 claim weak detections 
of Li with EWs between 30 and 70 mA in stars 6, 16, 20, 
59 and 126. Our measurements are in agreement apart from 



for star 6 (= JCO6_063) where we find EW[Li]= 99 ±31 mA 
compared with the 52 ± 5 mA measured by M08. 



4 LITHIUM DEPLETION IN IC 4665 
4.1 Lithium Equivalent Widths 

Figure [6] compares EW[Li] as a function of intrinsic V — I 
in IC 4665 (from Table 2 excluding the 12 non-members, so 
only a few contaminants may be left and these may have 
already been identified as the "possible" members - see sec- 
tion |3]) with EW[Li] measured in a similar way for stars 
in the IC 2391, IC 2602 and Pleiades clusters. The contri- 
bution of the weak, blended Fe I 6707. 44A line has not been 
removed from any of the EWs at this stage. We have initially 
chosen to make the comparison in this observational plane 
because the relative uncertainties are easier to understand 
and besides which, as the clusters have similar metallicity 
(see section [5l the transformation to Li abundance versus 
Tcff is similar for each of these clusters. 

IC 2391 has a measured LDB age of 50±5 Myr (Barrado 
y Navascues et al. 2004) and IC 2602 is usually assumed to 
have a similar age based on low-mass isochronal fits. The 
Li data come from the analyses of Randich et al. (1997) 
and Randich et al. (2001), the photometry from papers by 
Patten & Simon (1996) and Prosser, Randich & Stauffer 
(1996). The colours of the plotted data have been corrected 
for reddening {E{V - I) = 0.01 and 0.04 for IC 2391 and 
IC 2602 respectively - Patten & Simon 1996) . Our purpose 
is to show that whilst the amount of Li-depletion seen among 
F- and G-type stars (0.4 < {V — I)o < 0.9) is quite similar 
in IC 4665 and IC 2391/2602 (actually, it is consistent with 
very little depletion at all, see below and section |4^ . the few 
K- and early M-type stars we have found in IC 4665 seem to 
have undergone as much or more Li-depletion, despite this 
cluster having a younger LDB age. Note that uncertainties in 
the relative cluster reddenings (dominated by the ±0.06 mag 
in E{V — I) for IC 4665) affect the comparison of F/G stars 
far more than for cooler stars. 

The status of the K- and M-stars becomes more signifi- 
cant when considered in the context of the Li depletion pre- 
dicted by PMS models. Figure |6] also shows loci of EW[Li] 
expected for undepleted lithium, along with isochrones at 
ages of 30, 50 and 100 Myr, determined using the Baraffe 
et al. (2002) evolutionary models with convective mixing 
length set to 1.0 pressure scale height. These isochrones were 
calculated using a relationship between {V—I)o and effective 
temperature from Kenyon & Hartmann (1995), a curve of 
growth for the Li I 6708A feature described in Jeffries et al. 
(2003) and an assumed initial Li abundance of A{Li)= 3.3, 
on the conventional logarithmic scale (Anders & Grevesse 
1989). The curve representing an undepleted level of Li is 
an excellent match to EW[Li] measurements for cool objects 
in star forming regions over the whole colour range consid- 
ered (e.g. see Jeffries et al. 2009). The lower EW[Li] in stars 
of the same colour in young PMS clusters like IC 2391/2602 
has long been cited as evidence for PMS Li depletion (e.g. 
Randich et al. 2001). 

Evolutionary models by different authors predict differ- 
ent levels of Li depletion as a function of time (e.g. Pin- 
sonneault, Kawaler & Demarque 1990; Pinsonneault 1997; 
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Figure 6. The EW of the Li I 6708A feature versus intrinsic V — I for the IC 4665 candidate members (filled symbols, LDB age 28 ±5 Myr) 
compared with those in (a) the IC 2391 and IC 2602 clusters (open symbols, with an LDB age of 50 ± 5 Myr) and (b) the Pleiades 
cluster (open symbols, with an LDB age of 125 ± 8 Myr). In both plots triangles indicate upper limits. Also shown are model isochrones 
of Li depletion derived from the Baraffe et al. (2002) evolutionary models. Despite having a significantly younger LDB age, stars with 
1 < {V — I)o < 2 in IC 4665 seem significantly more depleted than the majority of objects in the older IC 2391/2602 clusters, more 
similar to the bulk of objects in the much older Pleiades and way below the predicted depletion levels for a 30 Myr old population. 



D'Antona & Mazzitelli 1997; Siess, Dufour & Forestini 2000; 
Piau & Turck-Chieze 2002) . We have chosen the BarafTe et 
al. (2002) models because they reasonably match the aver- 
age Li depletion seen in IC 2391/2602 at their LDB age, 
although underpredict the Li-depletion seen in some of the 
cooler stars. However, if IC 4665 were ~ 20 Myr younger 
than IC 2391/2602, as suggested by its LDB age, we would 
expect to see stars with much higher EW[Li] (at the same 
intrinsic colour) than in IC 2391/2602 for {V - I)o > 1, and 
to have detected lithium in stars as cool as {V — I)o — 2. 
Whether or not the detailed predictions of these models are 
correct, all published models predict a significant diminu- 
ition of EW[Li] with age at the colour of K- and M-stars. 
The observed progression of EW[Li] between the K- and M- 
stars (1.0 < (V - /)o < 2.0) of IC 4665 and IC 2391/2602 
is, if anything, the wrong way around. 

Figure [SJa shows a similar comparison with Li data 
from the Pleiades, which has an LDB age of 125 ± 8 Myr 
(Stauffer et al. 1998). The Li EWs come from the uni- 
form re-analysis by Sestito et al. (2005); upper limits for 
cool stars were added from Garcfa-Lopez, Rebolo & Martm 
(1994) and Jones et al. (1996). Kron V — I photometry from 
Stauffer (1982, 1984) and Prosser, Stauffer & Kraft (1991) 
has been converted to the Johnson- Cousins system using 
the transformation in Bessell & Weis (1987). A reddening 
of E{V — I) — 0.05 was assumed, except for those few 
stars listed in Table 9 of Soderblom et al. (1993a), which 
have larger reddening values. The F- and G-stars of the 
Pleiades have marginally lower EW[Li] than their equiva- 
lents in IC 4665. This may indicate some Li-depletion, but 
the uncertainty in E{V — I) for IC 4665 could eliminate 
most of this difference. The EW[Li] of the K- and M-stars 
of IC 4665 appear consistent with the most Li-depleted K- 
and M-stars in the Pleiades. Despite being much older, a 
significant fraction of the Pleiades population lies above the 
existing data in IC 4665. 

We have checked through all of our IC 4665 spectra 
to see whether any possible cluster members with high lev- 



els of Li might have been excluded during the RV selec- 
tion process. There are a number of objects with significant 
Li I 6708A detections that are either young field stars or 
candidate short-period binary stars in IC 4665 (see Fig. [3] 
and Table 1). However, none of these lie above the envelope 
of IC 4665 candidate members in Fig. |6l so their inclusion 
could only reinforce the phenomena described above. 



4.2 Lithium Abundances 

An alternative way to view these comparisons is to calcu- 
late Li abundances. The procedure described in detail by 
Jeffries et al. (2003) was used to convert intrinsic colours 
to Tcft and then transform EW[Li] into an abundance using 
a curve of growth for the Li I 6708A feature. The curve of 
growth is based on the calculations presented by Soderblom 
et al. (1993b) for Tcff > 4500 K and those of Zapatero-Osorio 
et al. (2002) for cooler objects. Prior to the transformation, 
EW[Li] was corrected for the presence of the weak, blended 
Fel line, using the formula EW[Fe]= 58 - O.OITch -I- 3.8 x 
10~^r^ff , derived from the correction as a function oi B — V 
suggested by Soderblom et al. (1993b) after applying the 
B - V/T^B relation of Kenyon & Hartmann (1995). The 
curves of growth are based on LTE atmospheres. NLTE cor- 
rections to the Li abundances were made using the results 
of Carlsson et al. (1994). The corrections vary from about 
-0.2 dex for Toff = 6000 K and A(Li)LTE = 3.4, to -f 0.25 dex 
for Tcff = 4500 K and ^(Li)LTE = 1.0. The Carlsson et al. 
(1994) corrections are only calculated for Teff > 4500 K, 
so for cooler stars no corrections have been applied. For 
Ted ^ 4000 K it is likely that the absolute corrections be- 
come smaller (Pavelenko & Magazzu 1996). 

The Li abundances for the IC 4665 members (and pos- 
sible members) from Table 2 are given in Tabled Note that 
any abundances A{Li)< —0.5 are extrapolations beyond the 
calibrated limits of the curves of growth and are pegged at 
this value. Similarly, we have not derived abundances for 
stars with Tcff > 6700 K. The abundance uncertainties are 
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Figure 7. A comparison of lithium abundances in IC 4665 (filled symbols) with those in (a) the IC 2391 and IC 2602 clusters (open 
symbols) and (b) the Pleiades cluster (open symbols). In both plots triangles indicate upper limits. Also shown are model isochrones of 
Li depletion from the Baraffe et al. (2002) evolutionary models at ages of 30, 50 and 100 Myr. The vertical bar at 3240 K indicates the 
approximate T^g of the LDB for IC 4665 found by M08. 



Table 3. Lithium abundances in IC 4665. Columns list the star, a V — I, the calculated T^g, the LTE Li abundance (j4(Li) = 12 + 
log[Af(Li)/Af(H)]) and its uncertainty and an NLTE Li abundance for stars with T^ff > 4500K (see text). 



Identifier 


V -I 


Teff (K) 


A(Li) 


o-A(Li) 


A{U) 


Identifier 


V -I 


Teff (K) 


A{U) 


o-A(Li) 


A(Li) 






(K) 


LTE 




NLTE 






(K) 


LTE 




NLTE 


JC01.427 


0.997 


5742 


3.33 


0.20 


3.12 


JCO5.280 


1.094 


5166 


< 1.82 




< 1.93 


JCO1.530 


0.787 


6389 


3.12 


0.19 


3.02 


JC05.282 


0.735 


6608 


3.51 


0.13 


3.34 


JC02.145 


0.768 


6450 


3.54 


0.10 


3.35 


JC05.296 


1.081 


5243 


3.17 


0.50 


2.96 


JC02_213 


0.933 


6062 


3.64 


0.12 


3.35 


JC05.329 


0.664 


7099 








JCO2_220 


2.023 


3856 


< 0.09 






JC05.394 


2.132 


3763 


< -0.37 






JC02.373 


0.883 


6214 


3.50 


0.09 


3.29 


JC05.472 


2.165 


3738 


< -0.50 






JC02.637 


2.745 


3456 


< 2.16 






JC05.515 


2.741 


3458 


< 0.56 






JC02.736 


2.986 


3393 


< 0.80 






JC05.521 


2.294 


3638 


< -0.50 






JCO3.065 


1.549 


4375 


0.75 


0.22 




JCO6.088 


1.028 


5587 


3.20 


0.17 


3.02 


JC03.285 


0.820 


6340 


< 3.33 




< 3.18 


JCO6.095 


1.013 


5668 


3.71 


0.35 


3.33 


JC03.357 


1.102 


5128 


3.20 


0.20 


2.97 


JC06.111 


1.780 


4120 


< 0.41 






JC03.395 


1.057 


5426 


3.41 


0.16 


3.13 


JCO6.240 


0.591 


7418 








JC03.396 


0.993 


5758 


3.29 


0.16 


3.09 


JCO7.021 


1.019 


5635 


3.37 


0.16 


3.13 


JC03.724 


2.310 


3625 


< 1.52 






JCO7.079 


0.847 


6286 


3.20 


0.13 


3.08 


JCO3.770 


2.699 


3470 


< 0.53 






JCO7.088 


1.882 


4001 


< 0.42 






JCO4.053 


1.164 


5071 


2.73 


0.15 


2.66 


JCO7.670 


2.576 


3489 


< -0.05 






JC04_226 


1.000 


5730 


3.36 


0.17 


3.13 


JC08.257 


1.051 


5459 


2.49 


0.19 


2.50 


JC04.337 


1.383 


4584 


< 1.06 




< 1.32 


JC08.364 


2.622 


3478 


< 0.11 






JC04.437 


2.079 


3806 


< -0.02 






JC08.395 


2.516 


3507 


< 0.64 






JC04.459 


2.577 


3488 


< 0.27 






JCO8.550 


2.839 


3430 


< 2.52 






JC04.591 


1.877 


4007 


< -0.50 






JCO9.120 


1.101 


5132 


< 1.97 




<1.86 


JC05.179 


1.240 


4900 


2.63 


0.37 


2.59 


JC09.281 


1.421 


4521 


1.29 


0.23 


1.55 



derived from propagated uncertainties in EW[Li] combined 
in quadrature with errors due to random Tetr uncertainties, 
assumed to be ±100 K. There are also a number of system- 
atic uncertainties to consider which are discussed at length 
by Jeffries et al. (2003). Most of these, such as the adopted 
temperature scale and curves of growth, affect stars in aU the 
clusters considered and do not aher comparisons between 
them. However, there are systematic uncertainties at the 
±0.2 dex level which will affect the comparison with models 
of Li depletion and the assumed initial Li abundance also has 
some uncertainty. We used an initial ^(Li)= 3.3, on the basis 
of the meteoritic Li abundance in the solar system (Anders 
& Grevesse 1989). It seems reasonable to assume that the 



initial abundance is similar for all the clusters considered, 
because they have very similar metallicities (see section |5]), 
but there could still be systematic errors of ~ 0.1 dex in the 
absolute model abundances. 

Figure [7] shows comparisons between IC 4665, 
IC 2391/2602, the Pleiades and the Baraffe et al. (2002) Li 
depletion models in the A{Li)-TeB plane. Because the trans- 
formations are similar for each cluster then the points made 
in the last subsection by comparing EW[Li] are equally valid 
here - namely, that the cool IC 4665 members are as or 
more depleted than stars of the same T^g in IC 2391/2602 
and occupy the same region as the most Li-depleted stars in 
the Pleiades. There is some evidence that the hotter stars 
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Figure 8. An absolute V vs intrinsic V — I diagram for members 
of the IC 4665, IC 2391, IC 2602 and Pleiades clusters. Also shown 
are two model isochrones at 30 and 100 Myr from Baraffe et al. 
(2002). The undertainty on the distance modulus of IC 4665 is 
indicated. The distance modulus uncertainties for the other clus- 
ters are small in comparison. IC 4665 is definitely younger than 
the Pleiades in this diagram but could be as old or younger than 
IC 2391/2602. 



(Teff > 5000 K) in IC 4665 have depleted a little less Li 
than their counterparts in the Pleiades, but again we cau- 
tion that this conclusion rests on the adopted reddening for 
IC 4665. A reduction in E{V — I) of 0.06 mag, would make 
the hotter stars cooler by about 100 K and reduce j4(Li) 
by about 0.15 dex. At lower temperatures {TeS < 5000 K) 
the assumed reddening has less influence; any uncertainty 
shifts the points almost parallel to the trend of decreasing 
abundance with decreasing temperature seen in the data and 
models, and hence does not greatly alter the comparison be- 
tween one cluster and another or the comparison with the 
models. 

The comparison with the models shows again that 
IC 4665 appears to have suffered much more depletion than 
expected. Other models with greater convective efficiency 
(e.g. the Baraffe models with a mixing length of 1.9 pressure 
scale heights) could address this, but these models would 
also predict far too much depletion to account for the ob- 
served Li abundances in IC 2391/2602 and the Pleiades. 

The Teff below which Li re-appears on the other side of 
the "Li-chasm" can be estimated from the data in M08. The 
LDB occurs at / — A" ~ 2.7, corresponding to Tctt — 3240 K 
from the colour- Toh relations of Kenyon & Hartmann (1995). 
This appears to be in perfect agreement with the predic- 
tions of the 30 Myr Baraffe et al. (2002) isochrone, though 
it should be noted that there are probably uncertainties of 
100-200 K in our estimated LDB T^g. 



5 DISCUSSION 

5.1 The age of IC 4665 

The LDB age of IC 4665 quoted by M08 is 28 ±5 Myr, where 
the bulk of the uncertainty is due to uncertainty in the dis- 



tance (370 ± 50 pc) to the cluster. The re-analysis of the 
Hipparcos data by van Leeuwen (2007) improves matters 
only a little, where a mean distance to IC 4665 of 35Gtll pc 
is found. The LDB age is found from the luminosity of the 
brightest very low-mass stars at temperatures below the "Li- 
chasm" (see Fig. 7) which still retain their initial lithium 
(see M08) . Although age spreads of up to 10 Myr have been 
suggested in young clusters (Palla & Stahler 2000; Jeffries 
2007), there is no strong evidence in M08 for Li-depleted 
cluster members below the LDB luminosity that might in- 
dicate the presence of a signiflcant fraction of older stars in 
the cluster. In any case, a 10 Myr age spread would not alter 
any of our subsequent discussion and conclusions. 

The LDB age has two key advantages over other age 
estimation techniques. First, although it is dependent on an 
assumed distance, it is much less sensitive to distance than 
some other methods, notably fitting isochrones to either the 
high- or low- mass stars (see below). Second, the physics of 
Li-depletion in fully convective low-mass stars is thought to 
be well enough understood that little room for systematic 
error remains. Indeed, any errors that persist must be due 
to effects that are currently not being debated by the vari- 
ous "flavours" of current low-mass evolutionary models, as 
they all give the same LDB ages to within a few per cent 
for clusters at ~ 30 Myr (see Burke et al. 2004; Jeffries & 
Oliveira 2005). For these reasons we attribute the greatest 
confidence to both the absolute and relative ages of clusters 
determined using the LDB technique. 

We can now ask whether the absolute ages and relative 
age ordering of clusters derived using other techniques con- 
cur with the LDB ages? For IC 4665, the age we would 
deduce from the low-mass isochrones in Fig. [S] (given a 
±0.25 mag distance modulus uncertainty) is in the range 
20-50 Myr, probably younger than the IC 2391 and IC 2602 
clusters (see also M08) and certainly younger than the 
Pleiades. These absolute isochronal ages could be in error 
so it helps to look at relative ages. Figure [8] plots absolute V 
versus intrinsic V — I diagram for these clusters, using the 
points plotted in Fig. |S] and using intrinsic distance moduli 
of 6.05 mag, 5.95 mag and 5.54 mag for IC 2391, IC 2602 and 
the Pleiades respectively (Stauffer et al. 1997). 

From Fig. |8] we can conclude that IC 4665 is younger 
than the Pleiades, but depending on the exact distance (the 
uncertainty is indicated in the diagram) it could be younger 
than IC 2391/2602 (and consistent with the LDB age) or 
of a similar age. Note that any uncertainty in the redden- 
ing tends to shift points parallel to the cluster sequence, so 
is relatively unimportant compared with the distance un- 
certainty. There is some indication from the upper part of 
the diagram, where the 30 and 100 Myr isochrones converge, 
that the assumed distance of 370 pc for IC 4665 may be too 
small. Better agreement with the shape of the isochrones 
would be achieved with a larger distance (say 400 pc) , which 
would move the IC 4665 towards a younger isochronal age, 
in accord with its LDB age and younger than IC 2391/2602. 

The ages determined from the nuclear turn-off in the 
Hertzsprung-Russell diagram do not contradict the age or- 
der of these clusters determined from their LDBs. IC 4665, 
IC 2391 and IC 2602 were assigned nuclear turn-off ages of 
36 Myr by Mermilliod (1981) who also assigned the Pleiades 
an older age of 78 Myr. These ages are dependent on uncer- 
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tain interior physics such as the amount of convective core 
overshoot, but the age sequence should be secure. 

Other indicators of youth such as the fraction of G-stars 
which are rapid rotators and the colour/spectral type at 
which Hq first appears in emission are less precise (e.g. Jef- 
fries, Totten & James 2000). In both these respects, IC 4665 
appears as young as IC 2391 and IC 2602, although the 
statistics are poor. In the range 0.6 < {V — I)o < 1-1, where 
significant spin down might be expected between ages of 30 
and 50Myr, there are 16/32 stars in IC 2391 and IC 2602 
(of those plotted in Fig. [Gj that have «sini > 20kms~^, 
where vsini measurements come from the same papers as 
the Li data. In the Pleiades this fraction has dropped to only 
6/56, where we have taken vsini from Queloz et al. (1998), 
Soderblom et al. (1993a) or Stauffer & Hartmann (1987), in 
that priority order. The equivalent fraction among IC 4665 
candidate members is 5/15, so is similar to IC 2391/2602, 
but inconsistent with the Pleiades (excluded with 95 per 
cent confidence by a two-tailed Fisher test). In IC 2391/2602 
Ha emission first occurs at {V — I)o — 1 and is ubiq- 
uitous for (V — I)o > 1.3 (see Stauffer et al. 1997; Jef- 
fries et al. 2000). In IC 4665 Ha emission is first seen in 
a star with {V - I)o = 0.84 ± 0.06 and 15/19 stars with 
{V-I)o > 1.2±0.06 exhibit Hq emission (see Fig.g}. Given 
the likelihood that a few contaminating dwarfs remain in the 
sample, the four cool candidates exhibiting Hq absorption 
are possibly non-members as discussed in section [3] 

We conclude that the isochronal age, the rotation of so- 
lar type stars and the chromospheric activity are consistent 
with the LDB age of IC 4665 and support its ranking within 
the age sequence of young clusters. 

5.2 Strong Li depletion in the K- and M-stars of 
IC 4665? 

Given the general agreement between the sequencing of the 
isochronal, nuclear turn-off and LDB ages, and the likeli- 
hood that there is little theoretical uncertainty in the latter, 
we are left with the task of explaining why, in IC 4665, the 
Li-depletion pattern seen among its K/M-stars would indi- 
cate an age that is significantly older. Even though only a 
small number of stars have been observed, one would esti- 
mate that IC 4665 is at least as old as the Pleiades from the 
evidence presented in Figs.[6jD and[7jD. 

One "uninteresting" possibility that can be discounted 
is that we have not observed any members of IC 4665 with 
1 < (V — /)o < 2 at aU, and that the smaU Li I EWs in 
Fig. El have been measured in contaminating field stars that 
randomly have the cluster RV. There are three strong ar- 
guments against this: (i) the predicted number of field star 
contaminants (16 ± 3 miims the 12 definite non-members, 
minus the 4 "possible members") remaining in the whole 
candidate sample is less than the number of cluster candi- 
dates in this colour range; (ii) six candidates in this colour 
range exhibit Ha emission, which is a strong indicator of 
youth and rarely seen in field stars with these colours; (iii) 
the presence of Li-rich mid/late M-dwarfs beyond the LDB 
in IC 4665 (see M08) conclusively demonstrates that the 
cluster mass function is not truncated below G-dwarfs. The 
limits of 1 < (V — /)o < 2 correspond approximately to 
1.0 > M/Mq > 0.7 (for a 28Myr Baraffe et al. 2002 
isochrone). de Wit et al. (2006) used a photometric survey 



and statistical subtraction of background contamination to 
estimate that the mass function of IC 4665, within 0.5 de- 
grees of the cluster centre, is dN/dM ~ 80 at these masses. 
Given that our survey has covered 87 per cent of "eligible" 
cluster members, with a further ~ 13 per cent lost through 
our RV selection (because only 87 per cent will be within 
1.5 sigma of the mean), we would expect to have seen ~ 20 
cluster members with 0.7 < M/Mq < 1.0, which is compa- 
rable to (and certainly not lower) than the 10 candidates we 
have actually found. 

More interesting physical explanations could centre 
around the chemical composition and rotation rates in 
IC 4665. Theoretically, PMS Li-depletion among stars with 
M > 0.5 Mq should be extremely sensitive to radiative opac- 
ity in the outer envelope and hence stellar metallicity. A 
higher metallicity should lead to deeper convection zones 
with hotter base temperatures and hence more extensive 
Li burning (see Pinsonneault 1997). For instance Piau & 
Turck-Chieze (2002) show that a mere 6 per cent increase in 
average metallicity could increase the PMS Li-depletion in 
solar mass stars by a factor of two. Likewise, the models of 
Siess, Dufour & Forestini (2000) predict that a 0.9 Mq star 
with a metallicity of Z — 0.02 will have depleted lithium 
by 1.2 dex at an age of 30 Myr, but if the metallicity is in- 
creased to Z — 0.03 the star has a similar Teff but has 
depleted lithium by more than 5 dex! 

At ~ 30 Myr a 0.9 star has {V - I)o ==; 1.3 and 
Tcff ~ 4300 K according to the Baraffe et al. (2002) mod- 
els we are using. Figures [6] and [7] show that in order to 
place the K-stars of IC 4665 above those of IC 2391/2602 
or the Pleiades would require about 1 dex less Li-depletion 
and about 2 dex less Li-depletion in order to agree with the 
predictions of the corresponding Baraffe et al. (2002) model 
isochrone. Perhaps this could be explained if IC 4665 were 
modestly more metal-rich than IC 2602/2391, the Pleiades 
and the metallicity adopted in the Baraffe et al. mod- 
els? The available evidence does not wholly support or ex- 
clude this. Shen et al. (2005) spectroscopically determined 
[Fe/H]= -0.03 ± 0.04 for F- and G-stars in IC 4665, with 
the abundances of other elements also similar to solar values. 
Platais et al. (2007) found [Fe/H]= -h0.06±0.06 for IC 2391, 
whilst the Pleiades has [Fe/H]= — 0.10±0.02 and solar abun- 
dance ratios (Wilden et al. 2002), both determined using 
similar techniques. 

A recent paper by D'Orazi & Randich (2009) obtains 
even more precise values of [Fe/H]= —0.01 ± 0.02 and 
0.00 ± 0.01 for IC 2391 and IC 2602 respectively, with the 
abundance ratios for other elements close to solar values 
too. It does not seem likely therefore that IC 4665 is more 
metal-rich than IC 2391 and IC 2602 unless there is some 
large discrepancy in the abundance analysis methodology 
of Shen et al. (2005) and D'Orazi & Randich (2009). Both 
analyses used Kurucz model atmospheres and derive abun- 
dances difi'erentially with respect to the Sun. However, Shen 
et al. (2005) fit their spectra with synthetic models to si- 
multaneously derive stellar parameters (including Tcff) and 
abundances, whilst D'Orazi & Randich derive Tcff from pho- 
tometric indices and then obtain abundances from a more 
traditional equivalent width analysis. A possible source of 
systematic error is therefore the adopted temperature scale, 
but Shen et al. show that their fitted Teff values are in excel- 
lent agreement with values determined from the relationship 
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Figure 9. A comparison of IC 4665 witii the Pleiades in the 
A{L\)-T^ff plane, where fast rotators {vsmi > 20kms~^) are 
given larger symbols. The plot also shows evolutionary tracks 
(from Baraffe et al. 2002) for stars of 0.7 and 0.%Mq to demon- 
strate that they have almost completed their Li-depletion by ages 
of 30 Myr and their evolution in this diagram is governed mainly 
by increasing T^g (see text). 



between T^s and [B — V)o proposed by Alonso, Arribas & 
Martinez-Roger (1996) using E{B - V) = 0.18 for IC 4665. 
This relationship is in turn very similar to the TcS /{B — V)o 
relationship used by D'Orazi & Randich (2009) (within 30 K 
at 5500 K), so no inconsistency is indicated. It would be a 
valuable exercise to analyse the IC 2391/2602 spectra using 
the same approach adopted by Shen et al. (2005) to rule out 
any remaining possibility of systematic error. 

In summary, IC 4665 could be of order 0.1 dex more 
metal-rich than the Pleiades - hence explaining the rela- 
tive Li abundances of the K- and M-stars of these clusters. 
However, on the basis of the published evidence, IC 4665 
seems unlikely to be more metal-rich than the IC 2391 and 
IC 2602 clusters by more than a few hundredths of a dex. 
If this were responsible for the excessive Li depletion of K- 
and M-stars in IC 4665 with respect to the predictions of 
PMS models, then the influence of metallicity must be very 
strong indeed. The lack of significant differences between 
the Li abundance trends seen in a number of Pleiades-age 
clusters argues against such a strong metallicity dependence 
during the PMS Li-depletion phase (e.g. Soderblom et al. 
1993b; Jeffries, James & Thurston 1998; Jeffries & James 
1999; Barrado y Navascues et al. 2001). 

Another possibility is a dependence of Li-depletion on 
rotation rate. Several authors have remarked that the wide 
spread in Li i EWs and consequently deduced Li abundances 
among the K-stars of young clusters may be correlated with 
rotation rate. This finding is based mainly on the Pleiades 
and clusters of similar age (e.g. Soderblom et al. 1993b; Jef- 
fries et al. 1998; Randich et al. 1998), where it is found 
that the fastest rotators tend to have the strongest Li I fea- 
tures and the lower boundary of the v4(Li) distribution is 
defined mainly by the slowest rotators. This may indicate 
that the fast rotators have preserved more of their initial 
Li, but could also be caused by problems with the model 
atmospheres, associated with rotationally-induced chromo- 
spheric activity or starspots (Stuik, Bruls & Rutten 1997; 
Ford, Jeffries & Smalley 2002). 



It is notable in this context that only 1 out of 10 IC 4665 
candidates with 1 < {V — I)o < 2 has vsini > 20kms~^. 
In contrast, the equivalent fractions for the IC 2391/2602 
and Pleiades samples in Fig. [6] are 10/25 and 18/44. Ac- 
cording to a two-tailed exact Fisher test we can marginally 
reject the null hypothesis that the fraction of rapid rotators 
is similar in IC 4665 to IC 2391/2602 or the Pleiades at 
confidence levels of 88 or 92 per cent respectively. Balanced 
against this, one might expect 0.7-1.0 Mq stars to spin-up 
as they contract between 28 Myr and 125 Myr, but even ne- 
glecting angular momentum loss this spin-up should only be 
a modest factor of ~ 1.2 (Baraffe et al. 2002). 

Figure |9] represents this information in graphical form, 
where the fast rotators in the Pleiades and IC 4665 are iden- 
tified by larger symbols. The Pleiades fast rotators appear to 
have depleted less Li than their slower rotating counterparts 
and these in turn have similar levels of Li depletion to the 
few K- and early M-type stars in IC 4665. Our view of the 
IC 4665 Li-depletion pattern might change if some rapidly 
rotating K-star members were found with higher levels of 
Li, but as discussed above, our data should include 75 per 
cent of cluster members in our survey area. If the RV se- 
lection criteria were relaxed so that 87 per cent of cluster 
members would have been included, there are still no Li- 
rich stars to be found. Even if a rapidly rotating, Li-rich 
member should emerge from the few remaining unstudied 
candidates, we would still have to explain why there are K- 
and M-type stars in IC 4665 that appear as Li-depleted as 
the most depleted stars of similar types in IC 2391, IC 2602 
and the Pleiades. 

Finally we comment on the status of the PMS models 
themselves. There may be a number of modifications or ad- 
ditions required in order to match emerging measurements 
of masses and radii for eclipsing PMS binaries (e.g. see Stas- 
sun et al. 2004). Much of the discussion centres around the 
efficiency of convection and whether it is simply less effec- 
tive than usually assumed (both on the PMS and for main- 
sequence binaries) or whether it may be suppressed by the 
presence of dynamo-generated fields at the interface between 
radiative and convective zones (e.g. Chabrier, Gallardo & 
Baraffe 2007). 

The lack of progression with age in Fig. [7] between 
IC 4665, IC 2391/2602 and the Pleiades (and other clus- 
ters at a similar age) offers an important constraint, but also 
presents an interesting dilemma when considered in this con- 
text. The problem is not how to prevent further Li depletion 
between ages of 30 and 125 Myr, but to explain why stars of 
a given Li abundance do not get significantly hotter over this 
period. This is demonstrated in Fig. (9] where evolutionary 
tracks for stars of 0.7 and 0.8 Mq are shown. In the Baraffe et 
al. (2002) models with a relatively low convection efficiency 
(a mixing length of 1.0 pressure scale height), Li depletion 
is almost complete by an age of 30 Myr at these masses. The 
change from the 30 Myr through to the 100 Myr isochrone is 
driven by the development of the radiative core, the corre- 
sponding change in the average polytropic index of the star 
and its traverse across the Hertzsprung-Russell diagram to 
hotter temperatures along the Henyey track (Henyey, Lele- 
vier & Levee 1955). To prevent, or minimise this evolution in 
Tas requires more efficient convection (or increased radiative 
opacities). This is clearly seen in current models with higher 
convective efficiency or higher metallicity (opacity), where 
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there is much less difference between the 30 and 100 Myr 
^(Li)-re£F isochrones (e.g. see Fig. 9 in Jeffries et al. 2003). 
But paradoxically, this leads to much greater Li depletion at 
a given age and T^s, so that such models cannot match the 
observed data, unless there were some compensatory effect, 
for example in the atmospheric modelling, which made the 
models much cooler (by several hundred Kelvin). 



6 SUMMARY 

In this paper we have used a radial velocity survey to se- 
lect candidate members of the young cluster IC 4665. The 
candidates have been filtered according to a number of cri- 
teria and the final list (Table 2) contains 40 highly probable 
members and a few possible members, of which we estimate 
only a few could still be contaminating foreground dwarfs. 
Table 2 should contain 75 per cent of cluster members with 
11.5 <V<18 over the central square degree of the cluster. 

The nuclear turn-off and isochronal ages of IC 4665 
are consistent with its lithium depletion boundary age of 
28 ± 5 Myr. However, the amount of lithum depletion seen 
among the K- and early-M stars of IC 4665 is greater than 
expected when compared with similar stars in other clus- 
ters of known age, and if judged in isolation, would indicate 
an age as old as the Pleiades at 125 Myr. Despite the small 
number of stars observed, this finding has important im- 
plications for the use of Li-doplotion as an age indicator. 
It is precisely in the regime of small number statistics for 
loose kinematic groups or even individual field stars that 
this technique has been most widely used in the literature 
to estimate ages (e.g. Zuckcrman & Song 2004; Barrado y 
Navascucs 2006; Mentuch et al. 2008). Our results suggest 
that Li-depletion among young late-type stars is a very crude 
age indicator indeed, with at least one other parameter be- 
sides age determining the rate and overall level of PMS Li 
depletion. 

A possible explanation could be differences in metal- 
licity and hence radiative opacities between different clus- 
ters/stellar groups, leading to variation in convection zone 
depth and PMS Li depletion. If so, then even mctallicity 
differences similar in size to their systematic uncertainties 
could be very important. Only a very careful, homogeneous 
study of Li-depletion in young clusters of similar age, but 
differing metallicities could confirm and calibrate this ef- 
fect. Alternatively, there are hints that the rotation rates of 
K- and M-stars in IC 4665 are slower than in comparison 
clusters and this may contribute to their apparent Li deple- 
tion. Demonstrating a causal effect will be difficult given the 
small number of K- and early M-type stars in IC 4665. 
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